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We analyze the large-rV behavior of SU (N) lattice gauge theories with adjoint fermions by study- 

i-C ing volume-reduced models, as pioneered by Eguchi and Kawai. We perform simulations on a 

single-site lattice for Nf = 1 and Nf = 2 Wilson Dirac fermions with values of N < 53. We show 
for both values of Nf that in the large-rV limit there is a finite region, containing both light and 

lf~) heavy fermions, of unbroken center symmetry where the theory exhibits volume independence. 

Using large-rV reduction we attempt to calculate physical quantities such as the string tension and 
meson masses. 
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1. Introduction 

The idea of volume independence of large-Af gauge theories dates back to the paper by Eguchi 
and Kawai [1]. They have shown that, provided certain conditions hold, in the large-N limit the 
pure gauge SU (N) theory defined on a single-site lattice 1 (the so-called Eguchi-Kawai model) has 
the same Wilson-loop expectation values as the theory defined on arbitrary (including infinite) 
volume. 

The only non-trivial assumption of their prescription is the lack of spontaneous breaking of 
the Zjv center symmetry. This assumption was quickly shown to be false [2]. The center symmetry 
in pure gauge theory is spontaneously broken below some critical lattice size L cr i t as can be seen 
both by lattice calculations [2, 3, 4] and in perturbation theory [2, 5]. 

Several ways to fix Eguchi-Kawai reduction were proposed over the years. Let us quickly 
sketch the three particularly interesting approaches: 

1. Use twisted boundary conditions (the so-called Twisted Eguchi-Kawai or TEK model) [6, 7]. 
While the original choice of twists has been shown not to work [8, 9, 10] a different choice 
appears to overcome the problems [11, 12]. 

2. Always work with L > L cr i t (b) [4]. This idea known as partial reduction or continuum reduc- 
tion 2 , allows one to make simulations directly in the pure gauge theory at the price of having 
to work with a box of finite physical size (thus with L cr i t — > oo in the continuum limit). 

3. Use adjoint fermions to stabilize the center symmetry. This idea, inspired by the orbifold 
large-Af equivalences [15, 16], was proposed in Ref. [17] and is the basis of our analysis 3 . 

In this paper we discuss the single-site lattice model with adjoint fermions - the Adjoint 
Eguchi-Kawai (AEK) model. It is shown in Ref. [17] that by adding Nf > 1/2 massless adjoint 
Dirac fermions obeying periodic boundary conditions, the center symmetry remains intact at the 
one-loop level. Perturbatively, fermions with masses larger than ~ 1 / (aN) were expected to break 
the center symmetry [20, 21]. However several lattice calculations [22, 23, 24] found that the cen- 
ter symmetry is likely to be intact for rather heavy fermions of mass ~ 1 /a, both for Nf = 1 and 
Nf = 2. In Ref. [25] we argued that this result holds in the large-N limit for Nf = 2. In this pa- 
per we present evidence that this is also most likely the case for Nf = 1. This result was given a 
semi-analytic understanding (using arguments going beyond perturbation theory) in Ref. [26]. 

Our results are of considerable phenomenological interest. Volume reduction in the region 
of heavy fermions, where the dynamics is governed mostly by gauge degrees of freedom, is a 
realization of the original idea of Eguchi and Kawai. Also, the Nf = 1 AEK model can be connected 
by a chain of large -N equivalences to the Corrigan-Ramond large-A^ limit of two-flavour QCD 
[17, 22]. The Nf = 2 AEK model, on the other hand, is supposed to behave similarly (e.g. with 
respect to the /3 -function) as St/ (2) gauge theory with 2 adjoint flavours which is most likely a 
conformal field theory (see Ref. [27] and references therein). 

'To be precise Eguchi and Kawai discussed the U (N) gauge theory that however coincides with SU (N) as N — > °°. 
2 See also Refs. [13, 14] for analysis in three dimensions. 

3 There is also a related idea of trace-deformed reduction [18, 19] that however becomes rather complex when 
reducing more than one dimension. 
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Motivated by the success of the reduction we attempt to extract large-distance quantities such 
as string tension and meson masses from our simulation. We discuss the N scaling of the effective 
volume and its impact on practicality of calculations in reduced models. 



2. Range of quark mass where volume reduction holds 

We analyze the single-site SU (N) lattice model with adjoint Wilson fermions by means of 
Monte Carlo simulations. The partition function is 

JT = Jd[U,y, V] e (5, «"»+E/£i Vj) , (2. 1) 

where 5 gaug e is the single-site equivalent of the Wilson plaquette action: 

W = 2Nb £ ReTrUf.UyU^Ul (2.2) 

(b = being the inverse 't Hooft coupling) and Dw is the Wilson Dirac operator on a single site 
(with periodic boundary conditions in every direction): 



D w = 1-k£ [(l-7.)^ dj + (l + 7M)C/r j 



where the hopping parameter K is related to the bare fermion mass by 



(2.3) 



m = — -4. (2.4) 
2k 

Note that the critical value K c (where the fermion mass vanishes) is shifted from 1/8 to a larger 
value due to additive renormalization and only becomes 1/8 in the continuum limit. 

Our numerical simulations are performed using the Hybrid Monte Carlo algorithm (equipped 
with the rational approximation in the Nf = 1 case). We have performed extensive scans of the 
K — b plane to find the phase structure of the model (see Ref. [25] for details on the Nf = 2 case) 
reaching N <53. Our main observable for the detection of breaking of the center symmetry are the 
"open loops": 

K„ = ^Tr U" 1 U'2 2 U" 3 U'J 4 , with n M = 0, ±1, ±2, . . . (2.5) 

where U~" = U^ n . Those results were supplemented by histograms of phases of link eigenvalues. 

The results for Nf = 1 and Nf = 2 are qualitatively similar. In both cases we find that, at fixed 
b, when decreasing the fermion mass from infinity (?c = 0) we go from a phase with completely 
broken center symmetry (as in the pure gauge model) to phases with partial breakings of Z^. When 
K exceeds a critical value Kf(b) (or, equivalently, when m is smaller than rrif) we enter the so-called 
"funnel" where the center symmetry is unbroken and volume reduction holds. 

An important question is whether the funnel remains of finite width in lattice units as ,/V — > 00. 
In Ref. [25] we performed an extrapolation to the large-A^ limit at b = 1 and found 

K f {b = 1, N = 00) =0.0655(5), (2.6) 
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(b) N f = 1 



Figure 1: Large-yV extrapolations of Kf. Note that the errors are rather conservative as the main goal was 
the exclusion of the "closed funnel" hypothesis. 



as shown in Fig. 1(a). Since we know that > 0.125, this implies that the edge of the funnel 
lies below K c even at N = °°. If we attempt to force a fit with Kf — > 0. 125 when N — > °°, we find 
# 2 /d.o.f. = 23, as shown in the figure. We have now performed a similar calculation in the Nf = 1 
case, obtaining 

K f (b = 1, N = ~) =0.0937(3), (2.7) 

as shown in Fig. 1(b). Thus the funnel is narrower at Nf = 1. The fit assuming Kf — > 0.125 at 
N = oo now has # 2 /d.o.f. = 2.9, so, while the closure of the funnel is not definitely excluded, we 
view this possibility as unlikely (particularly in light of the fact that K c > 0.125 due to additive 
renormalization, making the fit poorer). 



3. Physical quantities: preliminary results & outlook 

In Ref. [25] we showed that the AEK model with heavy quarks reproduces the large-volume 
pure gauge values of the plaquette very well. The natural next step using volume reduction is to 
calculate some long-range quantity such as the string tension. This can be done by calculating 
the static quark potential from reduced Wilson loops (wrapped many times around the single-site 
lattice). However, at our values of N we found it impossible to calculate the quark potential at 
separations large enough to give a reliable estimate. This was primarily caused due to large l/N 
corrections, and was not an issue of statistical errors. 

An example of the observed unphysical behaviour caused by the l/N corrections is the slow 
rise of Wilson loops at large separations. In Ref. [25] we showed that this can be understood 
qualitatively by inspecting the strong coupling expansion. On the single site every link in the 
Wilson loop comes in a pair with its hermitian conjugate and thus the zeroth order contribution in 
the strong-coupling expansion does not vanish. 

One way to tackle this problem is to use 2 4 lattices. The zeroth order does disappear for 
the odd-sized loops. In Fig. 2 we see that indeed the size of finite-yV corrections is substantially 
improved. However, our values of N to date (N < 15) are too small to allow extraction of the string 
tension. 
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(a) All Wilson loops (b) Only odd-sized Wilson loops 

Figure 2: Wilson loops on 2 4 lattice. N = 10, b = 0.35, K = 0.1. 



Another calculation that one can perform on a single site is finding meson masses using the 
"quenched momentum prescription" [28]. Our first attempt used fundamental Wilson and over- 
lap valence quarks over the sea of heavy adjoint Wilson fermions - this can be thought of as a 
quenched approximation as we expect the heavy adjoint fermions to have a negligible influence on 
the propagators. The method of quenched momentum prescription is based on inserting additional 
U(l) factors into the temporal gauge links, different for the quark and antiquark. These can be 
interpreted as momenta. For the pion the propagator one has 

Jt x (p , m) = Tr { 7s D 1 (U 4 e ip ' 2 , m) y 5 D" 1 (U 4 e^ 2 , m) } . (3.1) 

One can then extract the mass of the pion directly in momentum space. Our first results show 
sensible behaviour of the pion mass as a function of K", however we also see rather large l/N 
corrections and the values of N we used are too small for a reliable large-N extrapolation. 

The beauty of reduction is that both the physical volume and the color degrees of freedom 
of the unreduced theory are packaged into the color degrees of freedom of the reduced theory. 
Thus it makes sense, at least qualitatively, to introduce an effective box size, L e s(N). (For further 
discussion of the interpretation of this quantity see Ref. [25].) An important practical issue for 
numerical simulations is how L e ff scales with N. The reasonable possibilities appear bracketed by 
two cases — L e ff ~ N l / A , as motivated by the orbifold equivalence, or L e ff ~ N x l 2 (as for the Twisted 
EK model)[29, 30]. (The most optimistic case L e g ~ N l seems to be excluded by the data [25]). In 
Ref. [25] we argued that the Dirac operator spectrum prefers L e g- ~ N l l 2 , i.e. the more optimistic 
scaling behaviour. The situation is, however, far from clear, with different observables giving 
different indications. We hope that comparisons with results from a 2 4 lattice may help elucidate 
the issue. 

Another matter is the size of the finite-A^ corrections on the observables. We clearly see 
&{\/N) corrections in the plaquette and the Polyakov loops [25] that are rather large. It it thus 
promising to use twisted boundary conditions (the so-called Twisted AEK model [26, 31]) which 
have only G{\ /N 2 ) corrections and automatically have L e g ~ N l l 2 scaling. 

Finally, it would also be very interesting to compare the meson masses obtained using quenched 
momentum prescription with a more standard calculation involving a single elongated direction. 
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Such a setup also allows calculating the glueball masses. 
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